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Abstract 

We theoretically investigate a possible idea to introduce magnetic impurities to a superfluid 
Fermi gas. In the presence of population imbalance (A'^ > N^^ where N^j is the number of Fermi 
atoms with pseudospin a =t)i)) we show that nonmagnetic potential scatterers embedded in the 
system are magnetized in the sense that some of excess t-spin atoms are localized around them. 
They destroy the superfluid order parameter around them, as in the case of magnetic impurity 
effect discussed in the superconductivity literature. This pair-breaking effect naturally leads to 
localized excited states below the superfluid excitation gap. To confirm our idea in a simply 
manner, we treat an attractive Fermi Hubbard model within the mean-fleld theory at T = 0. 
We self-consistently determine superfluid properties around a nonmagnetic impurity, such as the 
superfluid order parameter, local population imbalance, as well as single-particle density of states, 
in the presence of population imbalance. Since the competition between superconductivity and 
magnetism is one of the most fundamental problems in condensed matter physics, our results 
would be useful for the study of this important issue in cold Fermi gases. 

PACS numbers: 03.75.Ss, 75.30.Hx, 03.75.-b 
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I. INTRODUCTION 



Magnetic impurity effects have been extensively discussed in the field of metallic super- 



conductivity. Within t 
Abrikosov and Gor'kov 



le Born approximation in terms of magnetic impurity scattering, 
l| (AG) showed that the superconducting phase is completely de- 
stroyed by magnetic impurities when the impurity concentration exceeds a critical value. 
Even below the critical impurity concentration, they showed that magnetic impurities re- 
markably affect superconducting properties, leading to the gapless superconductivity 
where the superconducting order parameter still exists but the BCS excitation gap is absent. 
We note that such a pair-breaking effect is absent in the case of nonmagnetic impurities, 
which is sometimes referred to as the Anderson's theoremlsl. 

Shiba extended the AG theory to include multi-scattering processes beyond the Born 
approximation [4 1 . He clarified that magnetic impurities induce bound states below the 
BCS excitation gap A. Using the numerical renormalization group theory, Shiba and co- 
workers further extended this theory to include the Kondo effect 5|, l6| . They showed that 
the transition from the Kondo singlet to the spin doublet of magnetic impurity occurs at 
Tk/A ~ 0.3, where Tk is the Kondo temperature. 



The recently realized superfiuid Fermi gases in ^°K[7| and ^Li|8|-|l0| have the unique 
property that one can tune the strength of a pairing interaction by adjusting the threshold 
energy of a Feshbach resonance. Using this advantage, we can now study superfiuid proper- 
ties from the weak-coupling BCS (Bardeen-Cooper-Schrieffer) regime to the strong-coupling 
BEG (Bose-Einstein condensation) regime in a unified manner ll|-|2l| . In addition, since 
the cold Fermi gas system is much cleaner and simpler than metallic superconductors, the 
former is expected as a useful quantum simulator to study various phenomena observed in 
the latter complicated systems, without being disturbed by extrinsic effects. Indeed, the 
pseudogap phenomenon, which has been extensively discussed in high-Tc cuprates 22|-|27| , 
has been recently observed in ^°K Fermi gases 28|, l29| . Since the latter system is dominated 
by strong-pairing fluctuations, we can now study how the pseudogap is induced by pairing 



fluctuations in a clear manner 



30 



32|. Thus, if one could study magnetic impurity effects 



in superfluid Fermi gases, it would be helpful for further understanding of the competition 
between superconductivity and magnetism, which is one of the most fundamental problems 
in condensed matter physics. 
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In this paper, we theoretically discuss an idea to realize magnetic impurities in a super- 
fluid Fermi gas. So far, nonmagnetic impurities have been realized in cold atom gases by 
using another species of atoms|33] and laser light {ImI]. One difficulty in realizing magnetic 
impurities in a superfluid Fermi gas is that spins {a =t, l) in this system are actually, not 
real spins, but 'pseudo'-spins describing two atomic hyperfine states. Thus, magnetic im- 
purities must be also 'pseudo '-magnetic objects acting on these atomic hyperfine states. In 
this paper, to realize such pseudo-magnetic scatterers, we use nonmagnetic impurities as 
seeds. Although the nonmagnetic impurities do not destroy superfluid state (Anderson's 
theorem), the superfluid order parameter may slightly decrease around the impurities, be- 
cause of the slight suppression of particle density, as schematically shown in FigHJ^a). In 
this case, when '["-spin atoms are added to the system, since they cannot form Cooper-pairs, 
they would behave like polarized magnetic impurities, and destroy the superfluid order pa- 
rameter around them. Thus, in order to minimize the loss of condensation energy by this 
pair-breaking effect, as shown in FigJTJ^b), excess ^-spin atoms are expected to be local- 
ized around the non-magnetic impurities, because the superfluid order parameter around 
the impurities has been already slightly suppressed before the ^-spin doping. The resulting 
nonmagnetic impurities accompanied by localized ^-spin atoms may be viewed as magnetic 
impurities. 

To conflrm this scenario, in this paper, we consider a two-dimensional attractive Fermi 
Hubbard model with a nonmagnetic impurity. Although this simple lattice model is different 
from the real three-dimensional continuum Fermi gas system, we emphasize that the presence 
of lattice, as well as the low-dimensionality, are not essential for the present problem. Within 
the framework of mean-fleld theory at T = 0, we calculate the spatial variation of the 
superfluid order parameter, as well as local population imbalance, around the impurity. In a 
polarized Fermi superfluid, we conflrm that the nonmagnetic impurity is really magnetized. 
We also show that the superfluid order parameter is damaged by this magnetized impurity, 
leading to bound states below the BCS excitation gap. 

We note that a different idea to realize magnetic impurities has been recently proposed in 



Ref. 



351], where an impurity with ^ is used (where a„ is a scattering length between an 



impurity and an atom with pseudospin a). We also note that, under the assumption of the 
presence of a magnetic impurity, properties of bound states have been examined in Ref. 



The outline of this paper is as follows: In Sec. II, we explain our formulation to calculate 
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FIG. 1: (Color online) Schematic picture of our idea to realize magnetic impurities. V{x) is a weak 
nonmagnetic impurity potential embedded in a superfluid Fermi gas. (a) In an unpolarized Fermi 
superfluid, although the nonmagnetic impurity does not destroy the superfluid state, the superfluid 
order parameter A(x) is slightly suppressed around the impurity, because of the decrease of particle 
density by the impurity potential V{x). (b) When an '|"-spin atom (solid circle with t) is added 
to the system, this excess atom is localized around the impurity, so as to minimize the loss of 
condensation energy by the pair-breaking effect. As a result, the impurity becomes magnetic, and 
the superfluid order parameter is remarkably suppressed around it. 

superfluid properties around a nonmagnetic impurity. In Sec. Ill, we numerically confirm 
our idea about the formation of magnetic impurities in a polarized Fermi gas. We discuss 
bound states induced by the magnetized impurity in Sec. IV. In Sec. V, we briefly examine 
effects of a trap potential, as well as finite temperatures. Throughout this paper, we set 
h = = 1, and take the lattice constant unity. 

II. FORMULATION 

We consider a two-component Fermi gas, described by pseudospin a =t, i- Assuming a 
two-dimensional square lattice, we put a nonmagnetic impurity at the center of the system. 
For simplicity, we first ignore effects of a trap potential, which will be separately examined 
in Sec.V. The model Hamiltonian is given by 

H = -t J2 + h.C -UJ2 + J2 [^(^ ~ l^'^] ^h'^^ (l) 

(i,j),a i 

where is the creation operator of a Fermi atom with pseudospin a at the i-th site. 
^j,cr = c|o-Cj,cr is the number operator at the i-th site, —t describes atomic hopping between 
nearest-neighbor sites, and the summation {i,j) is taken over nearest-neighbor pairs. —U 
is an on-site pairing interaction between '|"-spin atom and 4,-spin atom. Since we consider a 



polarized Fermi superfluid, the chemical potential Ha depends on a =t, J,. The nonmagnetic 
impurity potential V{i) is assumed to have the form 

VC) V (2) 

{Rx — RxoY + {Ry — RyoY + 

where {Rx-, Ry) is the spatial position of the z-th lattice site. {Rxo-, Ryo) is the center of the 
impurity potential. 

We treat Eq. ([T]) within the mean-field theory at T = 0. Introducing the superfluid order 
parameter A{i) = f/(cj_|Cj^-f-) (which is taken to be real in this paper), as well as the particle 
density no-(i) = (iT-i.o-), we obtain the mean-field Hamiltonian -^mf as 



^MF = -t E KA- + ^-C-] -E^( 



i i i 

In Eq. ([3]), we have ignored unimportant constant terms. (We also ignore them in the 
following discussions.) Sz,i = fii^^ — rii^^ is the spin density operator, and 

Szii) =n^{i) -n^ii) (4) 

represents the local magnetization. When this quantity becomes finite around the impurity, 
the last term in Eq. ^ works as an Ising-type magnetic impurity scatterer. The effec- 
tive impurity potential V{i) = V{i) + VHartree(0 involves the inhomogeneous Hartree term, 
VHartree(^) = ''{U /2)[n^{i) — n|(i) — uq], whcrc no is the particle density far away from the 
impurity, /i = yUav + {U/2)no is an effective chemical potential, where /Xav = [/^t + Z^^l/^- The 
difference of the spin-dependent chemical potentials 

^ = ^K-mJ (5) 

works as an effective magnetic field in Eq. (|3]). 

Since Eq. ([3]) has a bi-linear form, we can write it in the form -^mf = ^''Hmf^- Here, 
$t = 

(^i,t' ^2,t' ' ' ' ' ^M,t' ^i.-l-' ^2,1, ■ ■ ■ 5 cm,i), where M is the total number of lattice sites. The 
Hamiltonian matrix Huf is chosen so as to reproduce Eq. ([3]). As usual, Eq. ([3]) can be 



5 



diagonalized by the Bogoliubov transformation, 



(6) 



72 



7m 
7m+i 
7m+2 

\^ 72M J 

The 2Af X 2M orthogonal matrix W is determined so as to diagonahze ^mf- The diagonal- 
ized mean-field Hamiltonian has the form 

2M 



HuF = II^i7i7j, 

3=1 



(7) 



where Ej is the j-th cigcn energy. 

The superfiuid order parameter A(i), as well as the particle density (nj^o-) are given by, 
respectively. 



2M 



A{{)^uJ2WijWM+i,je{-Ej), 



2M 

3=1 

2M 

3=1 



(8) 
(9) 
(10) 



where Q{x) is the step function. The chemical potential /i^ is determined from the equation 
for the total number of cr-spin atoms. 



(11) 



To examine single-particle properties around the impurity, we consider the local density of 
states (LDOS). To calculate this, we introduce the single-particle thermal Green's function 
for 'I'-spin atoms, 

z, lUn) ^- dTe^'^^(r,{Q,t(T)cJ_^(0)}), (12) 
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FIG. 2: (Color online) Superfluid density of states in an unpolarized uniform BCS superfluid. 
The solid line and dashed line show the t-spin component p^{uj) and ^-spin component PjJ(w), 
respectively. In this figure, we have use the same values of the superfluid order parameter A and 
chemical potential fj, as A(l, 1) and fi^ in FigHJ respectively. 

as well as the Green's function for ^-spin holes, 

Gl{z,z,zun) = -jyTe^-^{TAcl{r)cM})- (13) 

Here, Un is the fermion Matsubara frequency, and /3 = 1/T is the inverse temperature. In 
the mean-field theory, Eq. f|T2l) is given by 

2M p 

G^^^ii,i,tun) = - W,,,W,,,, / rfre*"^(T.{7,(r)7j,(0)}) 

iy.2 

= E— (14) 

Executing the analytic continuation in Eq. ( !T4|) . we obtain LDOS for the ^-spin component 

as 

p^Uoo) = Im[GL(i, i, icOn u + i6)] 



2M 

EWt/(a;-E,). (15) 



In the same manner, LDOS for the J,-spin component is obtained from the analytic contin- 
uation of Eq. (|T3l) , as 

pU-) = -^Im[Gi,(.,M-„^c. + .5)] 
7 



2M 

Y,Wl_,,^^5{uj-E,). (16) 



In the absence of nonmagnetic impurity and population imbalance, Eq. (1151) reduces to 
the ordinary superfluid density of states in the uniform BCS state, 

P^^i^) = E^Siu - Ep) + + ^p), (17) 

p p 

and Eq. ( lT6l) also reduces to the BCS expression, 

pK^) = E«p^(^ + ^p) + E^p'^(^ - ^p), (18) 
p p 

where 



2 1 

U_ = - 

P 2 

2 1 



^^1 

-1-i^ 
- ^P 



(19) 
(20) 



In Eqs.( fT9l) and ( l20i) . i?p = + is the Bogoliubov excitation spectrum, where A is the 
uniform superfluid order parameter. = Ep—fi is the kinetic energy ^p = — 2t[cosp2;+cosp2/] 
of the tight- binding model, measured from the chemical potential /i. As shown in Figj2l ^^(0;) 
and p^{u}) have the BCS coherence peaks at a; ~ ±1.8t, as well as the van Hove singularities 
associated with the two-dimensional square lattice model at w ^ ±3.2t. From Eqs. ( !T7|) and 
f lTSj) . one finds that p^^{uj) describes '['-spin particle (hole) excitations when a; > (w < 0). 
p^l(a;) describes J,-spin hole (particle) excitations when u > {u < 0). 

We note that these physical meanings of p^-^{u}) and p^|(w) are slightly altered in the po- 
larized case {fif 7^ fii). In this case, within the neglect of the phase separation phenomenon, 
we again obtain Eqs. (fT7j) and (ITS!) . where the energy u is replaced hj co + h (where h is 
given in Eq. ([5])). Thus, in the polarized case, we need to replace the energy 00 hj u + h in 
the above discussions about the physical meanings of pf-|^(w) and p^j^(w). 

We also note that Eqs. ( |T5l) and ( fT6i) indicate that W^j and respectively, describe 

the probability densities of '|"-spin and J,-spin components of the j-th eigen state at the 
2-th site. When we write the wavefunction of the j-th eigenstate in the form \E'j(i) = 
(\&^(i), \I'^(i)), where each component is related to W^j and W^^^^j as 




FIG. 3: (Color online) Model nonmagnetic impurity potential V(i) used in this paper. We take 
{Ro^,Roy) = (21,21), Vo/t = 0.25, and T = 1. 

this wavefunction is normahzed as (Note that W is an orthogonal matrix.) 

M 2M 

E[i^?(or+i^?«r] = Ewt, = i. (22) 

j=l i=l 

From the physical meanings of p^-^{u}) and p^|(u;), one may interpret "^^{u)) and \l'^(z) as the 
t-spin particle and J,-spin hole components, respectively, when the eigen energy u is larger 
than —h. When u < —h, "^^{u)) and \l/j(i) have the physical meanings of '|"-spin hole and 
J,-spin particle components, respectively. We will use these interpretations in Sec. IV, when 
we examine physical properties of bound states. 

Before ending this section, we summarize details of our numerical calculations. We con- 
sider a 41 X 41 square lattice, and impose the periodic boundary condition. For the magnitude 
of on-site pairing interaction U, we take U/t = Q. To avoid band effects originating from 
the nested Fermi surface near the half-filling, we consider the low density region, by setting 

= = 300 in the absence of population imbalance (which gives the the particle density 
per lattice site [N^ + Ni]/M = 0.357 ^ 1). Starting from this unpolarized case, we gradually 
increase the number of ^-spin atoms to explore the possibility of localization of excess atoms 
around the nonmagnetic impurity at {Rxo,Ryo) = (21,21). For the impurity potential, we 
take Vo/t = 0.25, and F = 1 in the unit of lattice constant. The resulting impurity potential 
V{i = {Rx, Ry)) is shown in FigJ31 

Under these conditions, we numerically diagonalize -^mf to obtain W in Eq. ([6]), for a 
given set of %(i), %(«), /if; f^i)- We update them by using Eqs. (l8|)- (fTT|) . until the 

self-consistency is achieved. We then calculate superfiuid properties around the impurity, 
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such as LDOS's in Eqs. (fT5|) and ( !T6|) . as well as wavefunctions of low-lying excited states 
in Eq. 



III. FORMATION OF MAGNETIC IMPURITY AND PAIR-BREAKING EF- 
FECT 

Figure S] shows the superfluid properties around the nonmagnetic impurity in the po- 
larized case {N^ = 301 > = 300). For comparison, we also show the results in the 
unpolarized case in Figj5] (A^^ = A^^ = 300). As expected, FigsJll^a)-(c) show that the non- 
magnetic impurity is magnetized in the sense that an excess '["-spin atom is localized around 
it. Then, the last term in Eq. (|3]) works as a magnetic impurity potential. Indeed, FigJU^d) 
shows that the superfluid order parameter A{i) is damaged around the impurity, as in the 
case of ordinary magnetic impurity effect in metallic superconductivity [l|, |^. 

In the unpolarized case shown in FigJSl such a strong depairing effect is not obtained. 
The slight decrease of A{i) around the impurity is only seen in Figj5]^b), reflecting the slight 
suppression of the particle density n^{i) shown in FigjSj^a). This confirms that the remarkable 
suppression of A{i) in FiglD^d) is attributed to, not the original impurity potential V{i), 
but the localized ^-spin atom. The role of the weak nonmagnetic impurity potential V{i) 
is only to slightly decrease the superfluid order parameter around it, so as to capture the 
excess atom. 

When we increase the number of excess ^-spin atoms, they cluster around the impurit y, s o 



that local magnetization Sz{i) around the impurity increases, as shown in Figsl6](a)-(c)|37|. 
In addition, in the case of two excess atoms shown in Fig|6t^b), when we put one more 
impurity in the system, we obtain two magnetic impurities, as shown in Fig|6]^d), where 
each impurity is accompanied by one excess '|"-spin atom. Thus, using these, one can tune 
the magnitude of magnetization of each impurity, as well as the concentration of magnetic 
impurities. 

Here, we comment on the similarity between the magnetization of nonmagnetic impurity 



in FigJH and the phase separation observed in polarized Fermi gases 38|, l39| . In the latter 
phenomenon, the system spatially separate into the unpolarized superfluid region and po- 
larized normal region, and the minority J,-spin component is almost excluded from the latter 



region [38|-|45|. Such exclusion of J,-spin component from the polarized region can be also seen 
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FIG. 4: (Color online) Calculated superfluid state around the nonmagnetic impurity at {Rx,Ry) = 
(21,21). We take the population imbalance as N-^ = 301 > A'^ = 300, which gives the spin- 
dependent chemical potentials, fif = — 2.26t and m = — 5.20t. Panels (a) and (b) show the particle 
densities of t-spin atoms and |-spin atoms, respectively, (c) Magnetization Sz{i) = n-^ii) — n^{i). 
(d) Superfluid order parameter A(i). 
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FIG. 5: (Color online) Superfluid state around the nonmagnetic impurity at {Rx,Ry) = (21,21), 
in the unpolarized case {N^ = = 300). In this case, we obtain /^f = = — 3.73t. (a) no-(i) 
(cj =t)i)- (b) A(i). We note that Sz{i) = in the unpolarized case. 

in the present magnetic impurity problem. Indeed, FigJTl^a) shows that, while the density 
profile n^{i) of f-spin atoms is not so different from that in the unpolarized case (except for 
the enhancement at the impurity site {Ry = 21) due to the localization of the excess '["-spin 
atom), J,-spin atoms are clearly pushed out from the region around the impurity {Ry ~ 21). 

Another similarity between the present magnetic impurity formation and the phase sepa- 
ration in polarized Fermi gases is the negative superfluid order parameter {A{i) < 0) around 
the impurity seen in FigJTl^b). In a polarized Fermi superfluid, the possibility of Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) state (which is characterized by the spatial oscillation 
of the superfluid order parameter) has been proposed near the boundary between the super- 



fluid region and polarized normal region 46|]. The oscillation of the FFLO order parameter is 
the direct consequence of the formation of Cooper pairs with a finite momentum, originating 
from the mismatch of Fermi surfaces between the '["-spin component and J,-spin component 
in the polarized region. Although the polarized region is spatially small in the case of FigJTl 
one can still expect an FFLO-like pairing in the magnetized region, which leads to the sign 
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FIG. 6: (Color online) Intensity of the local magnetization Sz{i) around the impurity put at the 
center of the system. We take = 300. In panel (d), two impurities are put at {Rx, Ry) = (H, 21) 
and (31, 21). The spatial structure of each impurity potential is the same as Eq. ([2]). 

change of the superfluid order parameter seen in FigJTl^b). 

IV. LOCALIZED EXCITED STATES AROUND MAGNETIZED IMPURITY 

Figure E] shows the local density of states (LDOS) Pi,t(i^) around the impurity. In the 
unpolarized case (panel (a)), the impurity remains nonmagnetic, so that single-particle ex- 
citations are almost unaffected by impurity scatterings. LDOS far away from the impurity 
(solid line in FiglH]) agrees with the ordinary BCS density of states p^(w) shown in FigJ2J 
Even near the impurity, the overall structure of LDOS is unchanged (See the dashed line in 
FiglS), except for a slightly smaller excitation gap, reflecting the suppression of A(z) around 
the impurity seen in Figj5]^b). 

The situation is quite different in the presence of population imbalance. In this case, 
FigJH^b) shows that, in addition to the simple energy shift by the effective magnetic field 
h = — = — 1.47t, LDOS is remarkably modified around the impurity. Since the bulk 
BCS coherence peaks are at a; = — 3.32t and 0.38t (See the dashed line in FigJTOlfa).). the 
four peaks inside the BCS gap (— 3.32t < u < 0.38t) around the impurity seen in FigjHl^b) 
are found to describe low-lying excited states induced by the magnetized impurity. Since 
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FIG. 7: (Color online) (a) Particle density n„{i) {a =t;i) at = 21, as a function of Ry. We take 
= 301 and = 300. In this figure, na{i,N^ = N^) is the particle density in the unpolarized 
case (A'^ = = 300). (b) Spatial variation of the superfluid order parameter A(i) at Rx = 21. 
A(i, = N^) is the result in the absence of population imbalance (A'^ = = 300). 

these peaks disappear as one goes away from the impurity, they are bound states localized 
around the impurity. 

These bound states also appear in the J,-spin component of LDOS p^^(a;), as shown in 
Fig JTOT b). The fact that all the peak positions inside the BCS excitation gap (— 3.32t < a; < 
0.38t) are the same between panels (a) and (b) means that they are composite excitations 
consisting of particle and hole components. To confirm this in a clear manner, as example, 
we show the wavefunctions "^j^i) = (z), ^I/j (0) of the bound states (3) and (4) in FigJTTl 

In regard to the composite character of localized excited states, we note that the creation 
operator 7p ^ of the Bogoliubov excitation in the ordinary uniform BCS state is given by the 
sum of particle and hole creation operators as 

7i,t = ^pci.t + %c-p,;, (23) 
where and C-p,; are the creation operators of t-spin particle and 4,-spin hole, respectively. 

14 




P^((0) [t'] 

(a) A^i.=300, Ni=300 
1.2 r 




.10 -8 -6 -4 -2 2 4 6 8 10 



(b)A^^=301,A^i=300 
1.2 r 




.10 .8 -6 -4 -2 2 4 6 8 10 
CO/? 



FIG. 8: (Color online) Superfluid local density of states (LDOS) p^^{uj) around the impurity. In this 
figure, results at Ry = 21 show LDOS at the impurity site, (a) Unpolarized case {N-^ = = 300). 
(b) Polarized case (A'^ = 301 > = 300), where one excess t-spin atom is localized around the 
impurity. In calculating LDOS, we have added the small imaginary part 6 = 0.05i to eigen energies, 
in order to smear out unphysical fine structures around the coherence peaks, coming from discrete 
energy levels in the finite system. We briefly note that the influence of discrete levels still remains 
as sharp peaks around the van Hove singularity. We also use this prescription in FigsIUl and 1101 

That is, this basic property of ordinary single-particle Bogoliubov excitations is found to 
also hold in the bound states. 

As discussed in Sec. II, pf^^{uj) and p^_^(c<;), respectively, describe ^-spin particle (hole) and 
4,-spin hole (particle) excitations when u > —h {u < —h). Using this, we find that the 
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FIG. 9: (Color online) Local density of states /of^(w) in the unpolarized case (A''-^ = = 300). 
Apart from the fine peak structure seen around the van Hove singularity (which is due to the 
discrete levels in the finite system), LDOS far away from the impurity (solid line) agrees with the 
uniform BCS density of states p^{i^) in FiglJl 

bound states (1) and (2) in FigHni existing below —h = — 1.47t, are composed of J,-spin 
particle and '|"-spin hole excitations. Since the energies of the bound states (3) and (4) are 
larger than —h = — 1.47t, they are found to be composite excitations of f-spin particle and 
J,-spin hole. 

To see the mechanism of the bound state formation, it is helpful to examine the model 
Hamiltonian in Eq. (3) within the theory of magnetic impurity effect developed in the field 
of superconductivity 1, [3, Q-O], where the spatial variation of the superconducting order 
parameter is usually ignored. In addition to this simplification, when we also ignore the 
unimportant nonmagnetic potential V{i), and approximate the magnetic impurity scattering 
to a 5-functional potential, Eq. ([3]) reduces to, in momentum space, 

^ = E[^p - Hci,-Cp,- - ^ E [cp.t^-P,; + ^-c-] + ^ci,-Cp,a- (24) 

PiO^ P P,P',cr 

Here, Sz is the number of excess ^-spin atoms localized around the impurity. When Sz can 
take both positive and negative values depending on the direction of impurity spin, Eq. 
( 12^ is just the ordinary Hamiltonian describing a classical spin in a superconductor 
Although the pseudospin in the present case always points to the +2;-direction (when A^-^ > 
A^^), this difference is not crucial for the bound state problem. Indeed, Eq. gives the 
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FIG. 10: (Color online) Local density of states in the polarized case {N^ = 301 > = 300). (a) 
p^^{uj). (b) p^^{uj). The BCS coherence peaks exist at a; = — 3.32t and 0.38t far away from the 
impurity. (See the dashed lines in panels (a) and (b).) 

same bound state solutions as those obtained in the ordinary magnetic impurity problem j4|, 

,_,±Ai^M2MZ2!. (25) 

l + M0)s.C//2|2' ^ ' 

where we have approximated the normal state density of states to the value p(0) at the Fermi 
level, for simplicity. (For the derivation of Eq. fl2^ . see the Appendix.) While four peaks 
exist inside the BCS gap in FigHOl Eq. ( 125|) only gives two bound states. This means that 
the spatial variation of the order parameter ignored in Eq. flM|) is crucial for the formation 
of bound state in the present case. 

The superfiuid order parameter A(z) is known to have a property similar to the ordinary 



potential [47|, |48|]. Thus, the well structure of A{i) shown in Fig.(jl])(d) makes us expect that 
the bound states seen in Fig|TD]are trapped by this circular potential well|4i3]. Indeed, when 
we examine the bound states with u > —h, while the wavefunction of the lowest state (3) 
in Fig{TO]is almost isotropic and non-degenerate (See Fig JTW a).). the second lowest states 
(4) have the p-wave symmetry and are doubly degenerate, as shown in FigsJT2](bl) and 
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FIG. 11: (Color online) Calculated intensity of the wavefunctions of bound states. The upper and 
lower panels show the t-spin particle and J,-spin hole components, respectively. Panels (a) and (b) 
show the bound state (3) in FigHOl Panels (c) and (d) show the bound state (4). We note that 
the right panels show one of the two degenerate bound states existing at the peak (4) in FigJlOl 
For the other bound state, see Fig(T2{b2). 

(b2), which is consistent with the well-known result that the angular dependences of bound 
states in a circular potential well are given by cos(n6' + a) or sin(n6' + a) (n = 0, 1, 2, ■ ■ ■). 
In Figs JT2f bl) and (b2), one finds a = 7r/4, reflecting the breakdown of the rotational 
symmetry by the background square lattice. 

We also obtain localized states with higher angular momenta, such as the cZ-wave {n = 2) 
and /-wave (n = 3) states, as shown in the lower four panels of Fig{T2] (although they do 
not appears as peaks in Fig|10[50i]). Since their energies are close to the threshold energy 
u = 0.38t of the continuum spectrum, their wavefunctions spatially spread out, compared 
with the s-wave bound state in panel (a). Figure [12] indicates that the magnetized impurity 
would be also useful for the study of a quantum dot in a superfluid Fermi gas. 
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FIG. 12: (Color online) Intensity of the bound state wavefunction with uj > —h. Panel 

(a) shows the wavefunction of the bound state (3) in FigllOl The peak (4) in FiglTOl consists of 
doubly degenerate p-wave states given by panels (bl) and (b2). We also show other bound states 
with higher energies in panels (c) and (d) (although they do not appear as peaks in Fig 
The reason why the degeneracy is lifted in panels (cl) and (c2) is the broken rotational symmetry 
by the background square lattice. 

V. EFFECTS OF TRAP POTENTIAL AND FINITE TEMPERATURE 

So far, we have ignored effects of a trap, as well as finite temperatures, for simplicity. In 
this section, we briefly examine how they affect the localization of excess atoms. 

Figure ITST a) shows the intensity of local magnetization Sz{i) in the case when the gas is 
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FIG. 13: (Color online) Intensity of local magnetization Sz{i) in a trapped superfluid Fermi gas 
with population imbalance. Panels (a) and (b) show the case of harmonic trap Vj.|,^p(i) in Eq. (|26p . 
The population imbalance is chosen as (a) {N-^, Ni^) = (150, 70) and (b) {N-^, Ni^) = (150, 80). Panel 
(c) shows the case of quartic trap V^^ti in Eq. ([27]) with {N^, N^) = (150, 80). Panel (d) shows the 
case of box-type trap Vtrap(«) in Eq. ([28]), where we set {N^, N^) = (251, 250) and A = 2. In all the 
cases, we take V^^^-p = 5t. For the impurity potential at the center of the system, we take Vo/t = 2 
and r = 1. 



trapped in the harmonic potential, 

V!S,(, = (fl., fi,)) = ) " + ) 1 ■ (26) 



Here, (-R^^o, Ryo) = (21, 21) is the center of the lattice, and -Rmax = 41 is the number of lattice 
sites in the x- and y-direction. In this panel, the finite intensity of Sz{i) can be seen around 
the impurity at {Rx,Ry) = (21,21), which means that the magnetic impurity formation 
discussed in the previous sections also occurs in a trap. 

However, when the population imbalance is lowered, all the excess atoms are local- 
ized around the edge of the trap, so that the impurity remains nonmagnetic, as shown 
in Fig JT5r b). Thus, the population imbalance must be large to some extent to obtain the 
magnetic impurity. 

Since the excess atoms should avoid the spatial region where the potential is very large, 
the magnetization of the impurity is expected to occur more easily, when a trap potential 
around the edge of the gas is steeper than the harmonic trap in Eq. (I26l) . Indeed, in the 
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FIG. 14: (Color online) Local magnetization at the impurity site Sz{i = (21,21)), as a function of 
temperature. We take N-^ = 301 and = 300. The inset shows the calculated superfluid order 
parameter A as a function of temperature in the mean- field theory. In calculating A in the inset, 
we set = = 300, and ignored impurity potential. 



case of FigIT3l(b). when we replace the harmonic potential V^^^p{i) by the quartic one [5lH53|. 



trap 



xO 



R„ 



_^ ^ Ry- R 



yo 



(27) 



RxO Rma.x RyO ' 

the localization of excess atoms is realized, as shown in Fig JTSf c). As a more extreme case, 
when we use the box- type trap [5 



T/box/ T/0 



trap ^ 



trap 



(28) 



excess atoms can be localized around the impurity even in the case of small population 
imbalance, as shown in Fig JTST d). 



Next we examine effects of finite temperatures within the mean-field theory 55| . Figure 
[m shows the local magnetization Sz{i) at the impurity site, as a function of temperature. 
We find that the magnetization of the impurity remains finite even at finite temperatures. 
However, as shown in the inset, the present value of the pairing interaction U = Qt gives 
the mean- field superfluid phase transition temperature Tc/ 1 ~ 1 in the absent of population 
imbalance. Thus, the sudden decrease of Sz{i = (21,21)) at T/t ~ 0.08 means that the 
temperature must be far below Tc to realize magnetized impurities. 

In FigJlU we see the hysteresis of magnetization, which is characteristic of the first order 
phase transition. Regarding the present magnetization phenomenon as a kind of phase 
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separation, we find that this hysteresis behavior is consistent with the phase diagram of a 
polarized Fermi gas Sg], Is^] , where the first order transition from the phase separating gas 
to the superfluid phase occurs far below Tc in the case of small population imbalance. 



VI. SUMMARY 

To summarize, we have investigated how to realize magnetic impurities in a superfluid 
Fermi gas. In a two-dimensional attractive Fermi Hubbard model, we have calculated spatial 
variations of the superfluid order parameter, particle density, as well as polarization, around 
a nonmagnetic impurity within the framework of the mean-field theory at T = 0. 

In the presence of population imbalance (A^^ > Ni), we showed that the nonmagnetic 
impurity is magnetized in the sense that excess ^-spin atoms are localized around it. This 
magnetized impurity behaves like an Ising-type magnetic scatterer, so that the superfluid 
order parameter is destroyed around it. This pair-breaking effect is similar to the magnetic 
impurity effect discussed in the superconductivity literature. 

As another similarity between the present pseudo-magnetic impurity and real magnetic 
impurity in a metallic superconductor, the both impurities induce low-lying excited states 
below the BCS excitation gap. However, while the bound state formation by the latter 
magnetic impurity is usually discussed within the neglect of the spatial variation of the 
superconducting order parameter, we showed that the local suppression of the superfluid 
order parameter around the impurity is important in the former case. Namely, bound 
states are trapped inside a circular potential well formed by the superfluid order parameter. 
Because of this, the wavefunction of each bound state behaves like cos(n^ + Q;) or sm{n6 + a) 
(n = 0, 1, 2 ■ ■■) as a function of angle 6. Thus, this magnetized impurity may be also viewed 
as a circular quantum dot in a superfluid Fermi gas. 

Since the competition between magnetism and superconductivity is one of the most fun- 
damental problems in condensed matter physics, our results would be useful for the study 
of this important topic by using superfluid Fermi gases. In this regard, we briefly note that 
the magnetic impurity discussed in this paper is a classical spin, because it does not have 
an exchange term (which is symbolically written as s+cr_ + s+cr_, where s± = s^; ± iSy and 
cr± = (Ta; ± iay represent an impurity spin and the spin of a Fermi atom, respectively). As 
a result, the Kondo effect is absent. Thus, to examine the competition between the Kondo 
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singlet and Cooper-pair singlet in cold Fermi gases, it is an interesting future problem to find 
out an idea to realize a quantum magnetic impurity with an exchange term in this system. 
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Appendix A: Derivation of the bound state energies in Eq. (I25p 



To discuss the magnetic impurity problem in a Fermi superfluid, it is convenient to 
introduce the 4x4 matrix thermal Green's function 



G{v,v\iun) = - I rfre'"^(T4^p(r)vl/;,(0)}), (Al) 

J 

where = (Cp |, cj, j^, c_p j^, c_p^) is the four- component Nambu field 2|,|4|. In the absence 
of magnetic impurity, Eq.f lAip has the form G{p,p' ,iuJn) = Sp^piGo{p,iUn), where 

1 



G'o(p,iw^ 



(A2) 



(iUn + h)- ^pP3 + ApiTs 

Here, pi and [i = 1, 2, 3) are Pauli matrices, acting on particle- hole space and spin space, 
respectively. 



The full Green's function in Eq. flAip is obtained by including all magnetic impurity 
scatterings, which gives ^ 



G{p,p',iujn) = Sp^p'Go{p.iujn) + Go{p,iujn)i{iujn)Go{p',iujr. 



(A3) 



Here, the t-matrix i{iujn) involves effects of magnetic impurity smatterings as 

t[lUn) = —Szr3+{^ — S^j T3F{lUJn) + [ — Sz) T^F^lUJn) H 

U 



1 - ^s^F(za;„ 



(A4) 



where F{iujn) = J2]iGo(k.,iujn)T3- We briefly note that, in the case of ordinary magnetic 
impurity problem where the impurity spin can take Sz = ±s, we need to take the average 
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over the direction of impurity spin, as (s^") = s^", and (s^""*"^) = (n = 0, 1, 2, ■ ■ ■). In this 
case, i{iujn) is given byj4| 



(A5) 



As usual, the factor F{iun) in Eq. (IA4I) can be conveniently evaluated by replacing the 
momentum summation by the energy integration. Assuming the particle-hole symmetry of 
the fermion band, and approximating the normal state density of states to the value p(0) at 
the Fermi level, one obtains {2, 4| 



F{zun) = -7rp(0) 



{i Un + h)T3 + Api 



(A6) 



Bound states are obtained as poles of the analytic continued t-matrix, i{iun — )■ u + i5). 
Using this, we obtain the equation for the bound state energies as 



= det 



2 •/A2- (a; + /i)2 



(A7) 



which gives Eq. ( 125|) . We briefly note that the bound state equation in the case of Eq. f lA5 
is given by 



= det 



(A8) 



which also gives Eq. 
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